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ABSTRACT 

We show that if the turbulent magnetic difTusivity used in solar dynamos is assumed to be 'quenched' 
by increasing toroidal fields, much larger amplitude and more concentrated toroidal fields can be 
induced by differential rotation from an assumed poloidal field than if there is no quenching. This 
amplification and concentration mechanism is weakened and bounded by j X B feedbacks on the 
differential rotation. Nevertheless, it is strong enough to contribute to the creation of ^ 100 kC 
toroidal fields near the base of the convection zone, perhaps in conjunction with the 'exploding flux 
tube' process. Such high fields are necessary for sunspots to occur in low solar latitudes. 
Subject headings: Sun: interior — rotation — magnetic field 
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1. INTRODUCTION 

There are substantial theoretical reasons for conclud- 
ing that sunspots and active regions arise when toroidal 
flux tubes in or near the solar tachocline erupt to 
the photos phere (Fan et al.' 119941: iSchiissler et al.Ml994t 
iCaligariet al.. .1995. 1998). Given the Sun's rotation 
rate, for these eruptions to occur in low solar lati- 
tudes where spots are found requires that the tubes in 
the tachocline have field st rength approaching 100 kC 
llChoudhuri fc GilmanllT98l . Such large fields are a po- 
tential problem for dynamo theory applied to the Sun, 
since they are locally 10 — 100 times the equipartition 
value when compared to the kinetic energy density of 
both convective turbulence and differential rotation at 
tachocline depths. 

In all interface and flux transport dynamos applied to 
the sun, the strong toroidal fields achieved at tachocline 
depths come from the shearing by differential rotation 
there of a much weaker poloidal field. These toroidal 
fields are of much higher amplitude there than in the 
convection zone above because the magnetic diffusivity 
is assumed to be much smaller than in the much more 
turbulent convection zone. Since this shearing is also 
likely to lead to a back reaction on the local differential 
rotation, thereby extracting energy from it to amplify 
the magnetic field, the peak field achieved by induction 
depends critically on the rate energy can be resupplied 
to the differential rotation by the convection zone above. 
To be effective, this energy must be resupplied in a time 
quite short compared to a sunspot cycle. Since helioseis- 
mic inversions for differential rotation in the neighbor- 
hood of the tachocline do not show significant changes 
in the differential rotation within a solar cycle, whatever 
resupply is needed must be taking place. Indeed, the 
relative constancy of the observed differential rotation 
supports the validity of kinematic type dynamo models 
for the solar cycle, but sets constraints on the efficiency of 
the resupply mechanism and/or the toroidal field ampli- 
tudes produced by the shear. The resupply mechanism is 
not understood in detail, and any additional mechanism 
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that helps produce concentrated high amplitude toroidal 
fields potentially reduces the amount of resupply needed 
and may allow higher toroidal fields to be produced for 
the same feedback. 

A promising mechanism for achieving super- 
equipartition field strengths that does not depend 
totally on differential rotation is that of rising flux 
tubes 'exploding' befo re they reach the photo sphere 
l|Moreno-Insertis et al.i ri995: Re mpel fc Schiissleil f2001). 
In this scenario, many explosions of tubes result in 
significant loss of mass from inside the tubes, so that the 
part of the toroidal field that remains in the tachocline 
becomes more concentrated. If this happens at the 
unseen beginning of a new sunspot cycle, then the 
concentration of flux can lead to sufficient field strength 
that subsequent tubes get to the photosphere before 
they explode, so they can emerge as active regions. 

Here we propose another, perhaps complementary, 
mechanism for concentrating toroidal magnetic flux, in- 
volving the effect of turbulence, and therefore turbu- 
lent diffusivity, being partially or even largely sup- 
pressed where the toroidal fields are strong. The ef- 
fects of this " ry-quenching" have been considered more 
broadly in dy namo theory previo u sly by , for example, 
iRiidiger et all l|1994D and iTobiasI l|199fi|l . but without 
our em phasis on locally strong amplifi cation of toroidal 
fields. iRiidiger fc Kitchatinovl l)2000|) included the ef- 
fect of " ?7- quenching" together with quenching of tur- 
bulent heat conductivity in a model for the decay of 
sunspots. We use a formulation quite similar to that 
commonly used in dynamo m odels for " a- quenching" 
l|DikDati fc Charbonneai] I1999D (We make no attempt 
to evaluate the relative merits of various formulations 
of the ?7-quenching with gui dance from MHD turbulen ce 
theory, such as is done in .Vainshtein fc Cattaneol l)1992|) ). 
Our results below show that this Ty-quenching can be an 
extremely powerful mechanism for concentration when 
one is solving only the induction equation, but can be 
damped considerably when feedbacks on the differential 
rotation are taken into account. Thus this mechanism 
is subject to similar energetic limitations as when no Ty- 
quenching is allowed. 

2. SIMPLE MODEL 
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Fig. 1. — Schematic view of the geometry of our domain. In- 
dicated are also the imposed the velocity field u and the poloidal 
magnetic field Bp. 



There are several levels of complexity that could be 
used to illustrate how 77-quenching works. We choose 
here to work primarily in Cartesian geometry, since that 
allows us to obtain exact analytical solutions to illustrate 
the effect. We have also used a Cartesian finite difference 
code to compute the time dependent solutions of the in- 
duction equation and the reaction of the induced j X B 
force on the differential rotation. 

2.1. Governing Equations 

For our model, we consider a channel that has an x — z 
cross section, infinite in the down channel coordinate. 
Sidewalls are at a; = 0, L, bottom and top at z = 0, 
H. Roughly speaking, the x coordinate corresponds to 
colatitude on the Sun, the y coordinate, longitude and 
the z coordinate the vertical. The "pole" would be at 
X = 0, the "equator" at x = L. The geometry of the 
system is sketched in Figure ^ 

For the form of ry-quenching, we assume 
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in which 770 is the unquenched value, B = B{x, z) is the 
toroidal field to be solved for, and Bg is a constant, the 
value at which we assume 77 has been reduced by 50%. 

We seek solutions for the toroidal field B{x,z,t). The 
induction equation for this problem is given by 



(2) 



in which u is the assumed differential rotation, and B. 
Bz is the assumed poloidal field. We require 

dB:^ dB 




dx 
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but allow Br^ and/or B^ to cross the "equatorial" and/or 
top boundaries. For simplicity we omit meridional circu- 
lation from this calculation. 

There is a key transformation of the diffusion terms in 
equation jSJ that allows us to find analytical solutions, at 



least in the case of an assumed steady state. In particular 
from Buringtons tables, page 38, 
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from which we can write 
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As a starting point, we seek steady solutions to equation 
If we define a = iari~^ {B / Bs) then equation jSJ 
reduces to 

92 92 
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a Poisson-type equation for specific B^, B^ and u 
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2.2. Assumptions for differential rotation an poloidal 

field 

There are many choices for poloidal field and differ- 
ential rotation. Trigonometric functions are particularly 
convenient and illustrative. We choose 



. = -L/osin('^)sin(^) 



Bx — Bp sin 
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in which Uo and Bp are constant. For these cases, u 
has a maximum at a; = i (the "equator") and is zero 
at a; = 2L/3 (30° "latitude") as well as at a; = (the 
"pole"), very similar to the linear velocity of rotation on 
the Sun. There is also no differential rotation with x at 
the bottom, and a maximum amount at the top (anal- 
ogous to the solar tachochne). The poloidal field B^ is 
a maximum at a; = 0, and zero at a; = L, correspond- 
ing to a dipole field crossing the equator. B^ is zero 
at X = 0, and a maximum at x = L. By substitution, 
Bx, Bz satisfy equation JSJ. All these fields are shown 
schematically in Figure ^ Obviously many other choices 
of poloidal field and differential rotation are possible, but 
these seem particularly relevant to the solar case. 

2.3. Solution 

Substitutions of equations Q into the right-hand side 
of © yields 
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in which i?,„ — U^L/rjo is a magnetic Reynolds num- 
ber. By inspection, equation ((SJ should yield separable 
solutions 

(9) 



a = g{x)sm{^^ , 



from which we get 

B ^ B, tan 



g{x) sin 



H 



(10) 
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The tangent function is particularly critical in determin- 
ing the form of the toroidal field B. The substitution of 
equation into equation ^ , and changing variables by 
letting X = 2L\/tt leads to 

d'.9 .(lV BpR^ 
dA^'H^J .9 = __[3smAcos3A 

— cos A sin 3A] (11) 

For the solar tachocline, we can argue that L/H ^ 1 
(typical values are around 10 — 30) to a first approxi- 
mation, so we can ignore the second derivative term in 
equation Hll|l (unless we need to include a boundary layer 
to satisfy a boimdary condition). In that case 




and therefore the full solution for B is 




Thus we have found exact analytic solutions to a nonlin- 
ear partial differential equation! 



2.4. Qualitative Interpretation 

The solution given in equation H13|l. a particular case of 
the more general solution given in equation (|10|) , reveals 
many qualitative features, evident even without plotting. 
Of primary interest is that the induced toroidal field is 
proportional to the tangent of a function, of latitude and 
height, the magnetic Reynolds number, the ratio of as- 
sumed poloidal field to the quenching toroidal field, and 
the 'aspect ratio' of the domain, which is small in the 
actual solar tachocline. As is well-known, when the ar- 
gument p of a tangent function is small, then tanp w p. 
In this limit, equation H13|l becomes independent of the 
//-quenching field Bg. But tanp approaches infinity as 
the argument p approaches 7r/2. So we have the pos- 
sibility of very large amplification of the toroidal field 
from an initially relatively small poloidal field, as well as 
compared to B^. 

From equation H13() the induced toroidal field always 
vanishes at the 'pole' at x = 0, as well as the 'equator ' 
at X ~ L, and the top and the bottom of the layer. This 
is true even though the poloidal field crosses both the 
top boundary at z = H and the side boundary at x = L. 
Therefore the toroidal field is contained within the do- 
main of induction, and peaks in amplitude at mid height, 
z = H/2. In addition, by inspection of equation ((T^ . we 
can deduce that the argument of the tangent always has 
a peak at x = 2L /3 the equivalent of 30° latitude on the 
Sun. Obviously this location is a function of our assump- 
tions, but it follows from a differential rotation similar to 



4 



Gilman & Rempel 



80 
60 

^ 40 

CD 

20 




t = 1 .0 years 



Bs = 1 kG 

Bs = 3 kG 

Bs = 4 kG 

Bs = inf 



10 20 
Z [Mm] 



30 








100 200 
X [Mm] 



300 



80 

60 

^ 40 

CD 

20 








t = 2.0 years 



/ ~ 



10 20 
z [Mm] 



30 



80 
60 

^ 40 

CD 

20 








t = 2.0 years 



\ 



100 200 
X [Mm] 



300 



80 



60 
^ 40 

CD 

20 



t = 4.0 years / \ 
/ \ 



; 



\ 



' / \ \ 



\ 



10 20 
z [Mm] 



30 



80 



60 
^ 40 

CD 

20 




100 200 
X [Mm] 



300 



Fig. 3. — Snapshots of time dependent solution after 1, 2, and 4 years. Shown are vertical cross sections at x = 200 Mm = 2/3 L 
(left panels) and horizontal cross sections (right panels) at z = 15 Mm = 1/2 H of the toroidal field. Each panel shows a solution with 
no »7-quenching (solid), 77-quenching with Bs = 4kG (dotted), ?7-quenching with Bs = 3kG (dashed), and 7y-quenching with Bs = 1 kG 
(dashed-dotted) . 



that of the Sun, and a 'dipole-hke' poloidal field, so its 
location near the latitude where spots are first seen in a 
new sunspot cycle may not be coincidence. Bs is always 
positive, so for positive poloidal field Bp, the induced 
toroidal field is positive everywhere in the domain, even 
though not all terms in the induction forcing function 
(the right hand side) of equation © contribute with the 
same sign. 

Given the above, then, we can expect larger amplifi- 
cation of the toroidal field in the interior of the domain, 
peaking at mid-depth and x = 2L/3 or about 30° lati- 
tude. There can easily be a finite thickness and finite x 
or latitude range where the argument of the tan exceeds 
7r/2. What does this mean? In the section we dis- 
cuss reaching the steady solutions represented by equa- 



tion (|f 3(1 via time-dependent simulations that start from 
a state with zero toroidal field. We argue here that such 
solutions on the 'other side' of 7r/2 are not attainable in 
a finite time; the domain of such 'unreachable' solutions 
would be defined by the contour in the interior of our do- 
main on which, for given Bp, Bs and Rm, the argument 
of the tan = tt/2. On this contour, in effect steady so- 
lutions do not exist, because they would take an infinite 
time to reach the infinite field value. 

But so long as the argument of the tan remains smaller 
than 7r/2 in the whole domain, steady solutions are well 
defined throughout. Even in this case, however, it is of 
interest to know as a function of various parameters how 
long it takes a time-dependent solution to 'spin up' to 
the steady state. Obviously solutions that take the order 
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of a sunspot cycle or longer to reach the steady state 
will generally not be realized in the sun, since the cycle 
progresses and the peak toroidal field moves in latitude. 

The differences in predicted toroidal field amplitude 
as functions of Bg, Bp, and Rm in equation H13|l are 
straightforward to understand. Since the larger is Rm, 
the lower is the unquenched diffusivity, so we should ex- 
pect larger toroidal fields to result from a given Bp. And 
obviously the larger is Bp initially, the larger will be the 
resulting toroidal field. But the toroidal field is larger 
the smaller is Bg, since Bs is in the denominator inside 
the tangent function and in the numerator outside. This 
is because for smaller Bs, the local diffusivity is reduced 
starting at lower toroidal field, and so locally the toroidal 
field is kept from diffusing away; hence it amplifies fur- 
ther. For low enough Bg, so strong is the induction that 
no steady state can be reached. 

2.5. Quantitative results 

Figure 121 shows solutions of equation (|f 3|l for the com- 
mon parameters Rm = 3000, H/L = O.f , and Bp — 1 kG. 
To illustrate the amplification and concentration effect 
of ?7-quenching we have shown solutions for B^ = oo, 6, 
5, 4.3, and 4 kG. The top panels show a vertical cross 
section at a; = 200 Mm (left) and a horizontal cross sec- 
tion in the middle of the domain (right). The bottom 
panels show contours of the toroidal field for the cases 
with Bs = oo (left) and Bs — 4.3 kG. The case with 
Bs = 4: kG shown in the top panels is close to the transi- 
tion toward unbounded solutions (the peak field strength 
is around 200 kG). Figure |21 shows clearly the significant 



amplification and concentration effect discussed qualita- 
tively in section [T^ before. 

2.6. Solar effects? 

For plausible values of all the parameters in equation 
(|13|) . will the amplification effect we have found be sig- 
nificant for the Sun, in particular the solar tachocline? 
In the tachocline, we estimate that 30ms~^ < Uq < 
100 ms~^; lO^cm^s^^ < ?7o < lO^^cm^s^^, from weak 
overshooting turbulence there; 3 x 10^ km < L < 10^ km. 
These values lead to 10^ < Rm < 10^. For a tachocline 
of thickness ^ 3 x lO^'km, 0.03 < H/L < 0.1. Plausi- 
ble ranges for Bp and Bs are 5 x 10^ G < i?p < 5 x 10^ G 
and 10^ G< Bs< 10"* G. Then 

4xl0-3<^|p (^:^W„<4xl03 (14) 
At: Bs \ L J 

so in the Sun there is a full range of amplification possi- 
ble, from extremely small to extremely large, depending 
on the values of the various parameters we choose. Given 
that the peak value of the trig functions inside of equa- 
tion H13(l is about 2.6, even the smallest argument of the 
tangent is ^ 10~^, with 7r/2 — 1.57 corresponding to in- 
finite amplification; mid-range values for all parameters 
in (|14ll leads to an argument of the tangent for the solar 
tachocline of ~ 10, a factor of 6 above the first infinity. 
Thus it is easy to have solar conditions for which the 
amplification is large. 

3. TIME DEPENDENT SOLUTIONS 
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Fig. 5. — Influence of feedback on differential rotation on amplification process. We sliow snapshots of the solutions after 4 years. In 
each panel the solid line indicates the solutions presented already in Figure|3] the dotted, dashed, dashed-dotted, and triple-dashed-dotted 
lines represent solutions including feedback through Lorentz force and a restoration time scale for the differential rotation of r = 1, 3, 6, 
12 months, respectively. 



There are at least two effects not included in the solu- 
tions in (|13|l that could significantly diminish the ampli- 
fication effect on time scales relevant to a sunspot cycle. 
One, that we have already mentioned, is that the steady 
state is not reached before the solar cycle toroidal field 
moves to a lower latitude. The second is that as the 
toroidal field grows, there should develop e. j X B reac- 
tion force which damps the differential rotation respon- 
sible for the induction. Here we explore both of these 
effects with time dependent numerical calculations. 

We first focus on time dependent solutions only solving 
the induction equation and neglecting the feedback of the 
magnetic field on the shear flow. To this end we solve 
equation (5) using a finite difference code. We tested 
our numerical setup by comparing the stationary solu- 
tion obtained through our code with the solution shown 
in equation (13). For this comparison we have chosen 
the parameters in our numerical setup such that the ar- 
gument of the tangent in equation (|13|) does not exceed 
7r/2 anywhere in the domain to ensure the existence of a 
stationary solution. 

Common parameters for the time dependent solutions 
discussed below are: L = 3 x 10^ km, H = 'Sx 10"^ km, 
rjo = lO^cm^s-i , Uo = 100ms-\ and Bp = 1 kG. 
With this choice of parameters we have Rm — 3000, 
H/L = 0.1 and from equation (13) it can be expected 
that the argument of the tangent exceeds tt/2 if Bg = 
4kG. Figure 01 shows snapshots of time dependent so- 
lutions after 1 year, 2 years, and 4 years. The panels 



on the left show vertical cross sections at x = 200 Mm, 
where the field peaks in the horizontal coordinate. The 
panels on the right show horizontal cross sections in the 
middle of the domain (z = 15 Mm). We show a solu- 
tions with no //-quenching (solid line) and solutions with 
a value of Bs of 4 kG (dotted), 3 kG (dashed), and 1 kG 
(dashed-dotted). We can achieve similar toroidal field 
amplification patterns with Bs set much higher simply 
by assuming a smaller 770 and therefore a larger Rm- As 
explained above the solution with Bs = 4kG marks the 
transition from asymptotically stationary to unbounded 
solutions. The solution with no 77-quenching reaches an 
asymptotic field strength of around 6 kG, whereas the 
solution with Bs — 1 kG reaches around 80 kG after 4 
years of amplification through the shear fiow. Even the 
solution with moderate quenching _Bs = 4 kG reaches af- 
ter 4 years a value more than three times larger than the 
reference solution with no ry-quenching. Therefore it can 
be expected that the influence of 7y-quenching becomes 
visible even on time scales of the solar cycle. 

As an additional effect, the 7y-quenching leads to a con- 
finement of the magnetic field around the peak value, 
leading to a magnetic layer of less extent in radius 
and latitude compared to the reference solution with no 
quenching. The effect of confinement is largest for val- 
ues of Bs close to the critical value for which only the 
most central part of the magnetic field gets significantly 
amplified {Bg = 3kG and 4kG), whereas strong quench- 
ing {Bs — 1 kG) leads to a significant amplification of a 
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Fig. 6. — Contours of differential rotation (left panels) and toroidal field (right panels). The dotted curves in the left panels are the 
poloidal field lines. The top panels show the reference solution with no feedback, the middle panels the solution with r = 3 months, and 
the bottom panels the solution with r = 12 months. All snapshots are after 4 years of amplification time. 



much broader field profile. 

Figure01shows contour plots of the toroidal field after 4 
years for the cases shown in Figure |21 The spacing of the 
contours is equidistant, ranging from to the maximum 
value of each particular case. The greater concentration 
of toroidal flux around the center of the layer is clearly 
evident for = 3 and 4 kG. 

The time dependent solutions discussed so far show 
clearly that a significant amplification of a 1 kG poloidal 
field to several tens of kG toroidal field can be achieved 
within a few years, if 77-quenching is taken into considera- 
tion. However, these solutions assume that the feedback 
on the differential rotation through the Lorentz force can 
be neglected. In order to address this feedback in time 
dependent simulations we solve together with the induc- 
tion equation an equation for the differential rotation of 



the form 



du u ^ uq 

at " ~ 



Attq \ ^ dx ^ dz 



(15) 



In equation H15|l we include a drag-type forcing term 
which replenishes the differential rotation from the con- 
vection zone above on a time scale r. uo is the differen- 
tial rotation in the case of no Lorentz force feedback. For 
the density g we use a value of 210 kg m^"^ to represent 
the overshoot region. Since there are currently no theo- 
ries for the differential rotation in the tachocline, which 
would help to determine a reasonable value for r, we treat 
r as a free parameter and study the dependence of the 
solution on particular choices of values in the range of 1 
month to 1 year. Di fferential rotation models for the en- 
tire convection zone (jMiesch et al.l2000t lB run fc Toomrd 
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120021 lRempe]|l2005j) typically predict a time scale of sev- 
eral years; however, this time scale should be shorter if 
only the tachocline is considered. 

Figure |31 shows the influence of the feedback on dif- 
ferential rotation for different values of r ranging from 
1 month to 1 year. In each panel is shown the vertical 
profile of the toroidal field at a; = 200 Mm; the solid line 
corresponds to the solutions shown in Figure 13 after 4 
years. Even if a very short time scale for r of 1 month 
is used, the peak field strength reached is reduced by a 
factor of 2 in the cases with quenching of 77 (the case 
with no ?7-quenching shows a much weaker dependence 
on T since the field strength is very weak anyway for this 
choice of parameters) . For the choice of r = 1 year the 
peak field strength barely exceeds 5 kG. However, the 
tendency that ry-quenching leads to stronger and more 
confined field is still visible in all cases shown. 

Figure shows the profile of differential rotation (left 
panels) and toroidal magnetic field (right panels) for the 
cases with Bg — 1 kG. The top panels show the reference 
solution with no feedback, the middle panels the solution 
with r = 3 month and the bottom panels the solution 
with T = 1 year. In the left panels solid lines indicate pos- 
itive values of u, dashed lines negative values. We over- 
plotted the field lines of the poloidal field to show how the 
contours of the velocity field get aligned with the poloidal 
field lines to minimize the induction effect (Ferraros's law 
of iso-rotation). As a consequence, the feedback through 
the Lorentz force moves the shear layer of the differential 
rotation upward close to the equator (x = 300 Mm) and 
downward close to the pole {x — Mm) , leading to a sig- 
nificant reduction of the latitudinal shear in the middle 
of the domain, where the toroidal field peaks. 

To summarize, ?7-quenching has a significant influence 
on the field amplification, even if the time available for 
amplification is limited to a few years to reflect the solar 
cycle variability. Including the feedback on differential 
rotation through Lorentz force leads to a significant de- 
crease in the peak field strength reached; however, the 
influence of yy-quenching remains visible in most cases. 
Unless time scales of less than 1 month for the replen- 
ishment of differential in the tachocline are assumed it 
seems impossible to amplify a magnetic field to 100 kG 
as inferred by simulations of rising flux tubes. Other 
amplification mechanisms that do not rely on differential 
rotation as energy sour ce, as for example the ex plosion 
of magnetic fiux tubes (Re mpel fc Sch iissler 2001) or the 
downward draining of their interiors, most probably also 
play an important role. 

4. CONCLUSIONS AND DISCUSSION 
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We have demonstrated with a simple Cartesian model 
that the damping of turbulent magnetic diffusivity by 
the growth of strong toroidal fields is a powerful mecha- 
nism for making such fields even stronger and more con- 
centrated. Even though the expected reaction of these 
stronger fields on the differential rotation that induced 
them can limit the amplitude of this effect, it should re- 
main an important contributor to the overall operation 
of the dynamo inside the Sun. 

Our result is qualitatively consistent with the effect 
common to all flux transport and interface dynamos for 
the sun, in which the high toroidal field in the tachocline 
is achieved because of the low magnetic diffusivity there. 
In that case, the low diffusivity comes from the assumed 
weak turbulence there, while in the present work the 
toroidal field is assumed to locally reduce t he turbulence . 
This quenching has already been shown bv lTobia^ 1)19961) 
to have profound effects on, for example, dynamo mode 
type, when applied to more idealized interface dynamos 
considered for the bottom of the convection zone. 

The role this toroidal field amplification mechanism 
plays in current dynamo models for the Sun should 
be tested. While our model used Cartesian geometry 
because of its simplicity and because it yielded useful 
analytical solutions, we have done sample calculations 
in spherical geometry that show similar effects will be 
present there. 

It will be interesting to see even in the kinematic dy- 
namo regime how much concentration and amplification 
of toroidal field occurs when a real dynamo solution is ad- 
vancing through a simulated solar cycle, so the toroidal 
field is continually moving toward the equator. This ef- 
fect could also be simulated with our simple Cartesian 
model by taking poloidal fields that represent various 
phases of a solar cycle; the case we have chosen corre- 
sponds roughly to the maximum phase, since the poloidal 
field peaks at the pole and there are no separate 'old' and 
'new' cycle poloidal fields present. It is clear from our 
basic equations that analytic solutions should also exist 
for these other cycle phases. An important factor limit- 
ing toroidal field amplification will be simply how long in 
reality it takes for the imposed poloidal field to move in 
X (latitude) enough to change significantly the location 
where induction of toroidal field is strongest. 
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